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I We examine the pion exchange mode of neutrinoless double beta decay 

{QuPP) induced by the i?-parity violating quark- lepton operators of the su- 
persymmetric (SUSY) extensions of the standard model of the electroweak 
interactions. The corresponding nuclear matrix elements are evaluated 
within the renormalized quasiparticle random phase approximation with 
^ ' proton-neutron pairing, which includes the Pauli effect of fermion pairs 

1/^ , and does not collapse for a physical value of the nuclear force strength. 

It is argued that the pion-exchange mode of Oi^/3/3-decay dominates over 
' the conventional two-nucleon mode in the case of the SUSY mechanism. 

As a result sensitivity of 0!//3/3-decay to the SUSY contribution turns out 
to be significantly better that previously expected from the two-nucleon 
mode calculations. An upper limit on the i?-parity violating coupling A'm 
is derived from non-observation of neutrinoless double beta decay. This 
' limit is much stronger than that expected from the near future accelerator 

experiments. 

Oh: 

' 1 Introduction 



The observation of neutrinoless nuclear double beta decay {Qvj3(3) would un- 
^ , doubtedly indicate the presence of new physics beyond the standard model 

' (SM) of electroweak interactions. However, as yet there is no any experimental 

evidence for this lepton-number violating exotic process. On the other hand 
non-observation of 0;//3/3-decay at given experimental sensitivity allows one to 
establish limits on some parameters of new physics. The most famous example 
is the light effective Majorana neutrino mass which is limited by the QvjSP-Aecay 
experiments [Q at the level (m'^j) < 0(0.5 — 1.1 eV) ||l|, || depending on nuclear 
model. 

It is known that the Majorana neutrino exchange is not the only possible 
mechanism of Oi^/3/?-decay. The lepton-number violating quark-lepton inter- 
actions of the R-parity non-conserving supersymmetric extensions of the SM 
{l^p SUSY) can also induce this process It is worthwhile to notice 

that the corresponding 1st generation I^p Yukawa coupling are so stringently 
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constrained by non-observation of the Oi^/?/3-process that all possible 1st gener- 
ation ^peffects are pushed beyond the reach of the present and the near future 
accelerator and the other non-accelerator experiments ^ . 

In searching for tiny effects of the physics beyond the SM the main dis- 
advantage of Oi^/3/9-decay experiments compared to the accelerator ones is the 
necessity of taking into account nuclear structure. 

Although there are many difficulties the study of Oi^/3/3-decay has the ad- 
vantage of unprecedented accuracy and precision with which the process can 
be studied experimentally. It is possible to observe large samples of several 
kilograms of potentially decaying nuclei and to search for the decay of a single 
nucleus. The lower half-life limit for Oi//3/3-decay measured in this way is very 
high. In addition, the 2i//3/3-deca.y process predicted by the standard model can 
be measured at the same time in order to check the nuclear structure calcula- 
tions. 

The nuclear Oi^/3/3-decay is triggered by the Qvj3[3 quark transition d-\-d u-\- 
u + 2e~ which is induced by certain fundamental interactions. It was a common 
practice to put the initial d-quarks separately inside the two initial neutrons of 
a Oj//3/3-decaying nucleus. This is the so called two-nucleon mode of the Qv(3l3- 
decay. If the above quark transition proceeds at short distances, as in 

the case of I^p SUSY interactions, then the basic nucleon transition amplitude 
n-\-n^p + p + 2e~ is strongly suppressed for relative distances smaller than 
the mean nucleon radius. 

The goal of this paper is to discuss the pion-exchange SUSY mechanism 
which is based on the one and two pion exchange between the decaying neu- 
trons. The two pion exchange counterpart of this mechanism was first studied 
in Ref. At the quark level this mechanism implies the same short-distance 
I^p MSSM interactions as in Ref. Q. However, it essentially differs from the pre- 
vious consideration of the SUSY contribution to the Oi^/3/3-decay at the stage 
of the hadronization. We will show that in the case of the I^p MSSM induced 
quark transition the pion-exchange contribution absolutely dominates over the 
conventional two-nucleon mode. As a result, a significant improvement of the 
previously known |^ two-nucleon mode Qv {3 (i-decay limit on -Yukawa cou- 
pling X'lii becomes possible. 

We calculate the nuclear matrix elements governing the two-nucleon, one 
pion-exchange and two-pion exchange SUSY contributions to Oi'/J/^-decay within 
the renormalized Quasiparticle Random Phase Approximation with proton- 
neutron pairing (fuU-RQRPA) . This nuclear structure method has been devel- 
oped from the proton-neutron QRPA method, which has been frequently used in 
the Oi^/3/3-decay calculations. The fuU-RQRPA is an extension of the pn-QRPA 
by considering the effect of the proton- neutron pairing and the Pauli effect of the 
fermion pairs in an approximate way. In this way the sensitivity of the nuclear 
matrix elements to the details of the nuclear Hamiltonian is reduced consider- 
ably and more reliable values on the lepton number non-conserving parameters 
are obtained. 
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2 The Nuclear Structure Method 



Most of the double beta decaying nuclei are not yet accessible for the detailed 
shell model treatment because of the complexity of this approach. Therefore the 
pn-QRPA is mostly used for the nuclear structure calculations. In general 
the half-life for a 0z^/3/?-process can be factorized in the form 

= G\ME\\ (1) 

J- 1/2 

with the leptonic phase space factor G, the nuclear matrix element ME and 
the factor e that is given by the special extension of the standard model under 
consideration. 

The nuclear matrix element ME is ruled by the nuclear structure of the 
involved nuclei, but also by the considered Oi//3/3-decay mechanism, which de- 
termines the transition operators. We calculate the matrix element in the in- 
termediate nucleus approach, which requires the construction of a complete set 
of the intermediate nuclear states. 

In the framework of the QRPA or renormalized QRPA the to*'' excited states 
with angular momentum J and projection M is created by a phonon-operator 
Q with the properties 

Q7LKpa)^\^^JM) and Q|0+p^)=0. (2) 

Here, |Oj^p^) is the ground state of the initial or the final nucleus. The phonon- 
operator Q takes the following form 

Q7L^ = E ^(M^).^^^(M^, JM) + y(^^,,),^i(M^^, JM) (3) 

where the tilde indicates time- reversal. A^^k^liy, JM) is the two quasi-particle 
creation and annihilation operator coupled to good angular momentum J with 
projection M namely 

A\k^ilv,JM)^n{k^iM E • (4) 

The indices /i and v denote the isospin structure of the bifermion operator. In 
the case of a BCS definition of the quasi-particles without pn-pairing it would 
distinguish between proton and neutron type quasi-particles. For the case of a 
HFB definition of the quasi-particles with pn-pairing the quasi-particles have no 
definite isospin anymore and the index becomes 1 or 2 ||]. The factor n(fc/i, Iv) 
guarantees the normalization for the case of two identical particles fc/i and Iv. 
From Eq. (|) and (|) one can derive the RQRPA equation 
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where fJjl denotes the excitation energy of the m*'' excited state with angu- 
lar momentum J and parity tt with respect to the ground state |Ojp^). The 
three matrices A, B, U are given by the expectation values of the following 
commutators in the correlated ground state 

Kt^J = K.pa\ [A{aahp, JM), [H,A\k^^l,., JM)]] |0+p^) , (6) 

= Kpa\ [A{aabp, JM), [A{kfil,., JM),h]] |0+p^) , (7) 

^fcT'i!^ = Kpa\ [Aiaabp, JM),A\k^ilv, JM)] |0+p^) . (8) 

For solving the above equation it is necessary to introduce some approxima- 
tion scheme for the calculation of the A, B and U matrices. The simplest and 
the most frequently used one is the quasiboson approximation (QBA) scheme, 
which implies the two-quasiparticle operator A^{kfiliy, JM) to be a boson op- 
erator. However, it is well-known that the QBA violates the Pauli exclusion 
principle because we have neglected terms coming from the commutator of the 
two bifermion operators by replacing the exact expression for this commuta- 
tor with the its expectation value in the uncorrelated BCS/HFB ground state, 
which is the vacuum for the quasi-particle operators. It turns out that the QBA 
is a poor approximation and leads to too strong ground state correlations, an 
unphysically large occupation of quasi-particle states in the correlated ground 
state, close to a collapse of the QRPA solution. But the ground state correla- 
tions influence the nuclear matrix elements severely and in general the use of 
the QBA leads to a very sensitive dependence of the nuclear matrix elements 
on the strength of the residual interaction in the particle-particle channel. 

To overcome this problem the Pauli-principle has to be incorporated in the 
approach by which the occupation of the quasi-particle states in the 

correlated ground state would be limited. The Pauli-principle can be incorpo- 
rated to large extend by calculating explicitly the commutator of Eq. (||) with 
the single quasi-particles obeying Fermi-anti-commutation rules. Neglecting the 
non-diagonal part the commutator is not anymore boson like, but obtains cor- 
rections to its bosonic behavior due to the fermionic constituents. 

K^il^j = {0+p^\[A{aabp,JM),A^ikf,li^,JM)]\0+p^) 

~ n{kn, lv)n{aa, b(3) (SkaSf_,aSibS^p - SiaSuaSkbS^,p{~l)^''^^'^'^ 

1 - ir(ojjPAl[«i:fc%fc]oo|o+p^) j' 



(9) 



with the known abbreviation j — ^ {2j + 1). This expression is still diagonal 
in the quasi-particle configuration indices, and therefore the QRPA Eq. (||) can 



4 



easily be brought to eigenvalue form 



As in this form the matrices A and S and the eigenvectors are renormalized by 
D this approach is called renormalized QRPA (RQRPA) . To solve this equation 
the quasi-particle occupation in the correlated ground state needs to be known. 
But the occupation can only be derived from the back-going amplitudes Y of 
the RQRPA diagonalization. An implicit formulation for the I?-matrix can be 
given, if one express the one-body densities of Eq. (^) in terms of mapping on 
two boson operators. Then one arrives at the formula 

D{kiii^)j = l - j^"^ D(^k,ik'^i')j'J''^\Y{ktj.k'tj.')j'^\ (11) 

J'm 
.J'm 

Eq. (|l^) and (|lO|) have to be solved self-consistently in an iterative proce- 
dure. Note that an eigenvalue equation has to be written down separately for 
every multipolarity J'" which are then coupled by Eq. (pi]). The method pre- 
sented above include both proton-neutron pairing and the Pauli effect of fermion 
pairs and is denoted the fuU-RQRPA. In the limit the proton-neutron pairing is 
switched off, one obtains the pn-RQRPA method. 

In order to calculate double beta transitions two fully independent RQRPA 
calculations are needed, one to describe the beta transition from the initial 
to the intermediate nucleus and another one for the beta transition from the 
intermediate to the final nucleus. For that purpose the one-body transition 
densities of the charge changing operator has to be evaluated. They take the 
following form: 

< J"^^ II [Cpfc^n;]./ II 0+ >= V2JTT X (12) 



< 0/ II [Cpk'~^ni'].j II J'^mj >= V2JTT x (13) 
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with m (//a, 1/5) = (14,^ g 'Ji "2'' ■ We note that the J'^) and J'^) 

amphtudes are calculated by the renormalized QRPA equation only for the 
configurations fj,a < i/b { i.e., fi = v and the orbitals are ordered a < b and 
fj, = 1, u = 2 and the orbitals are not ordered). For different configurations 
X^Jk, I, J"") and F"^(fc, I, J'^) in Eqs. (ph and (Q are given by following the 
prescription in Eqs. (65) and (66) of Ref . [pi] . The index i (f) indicates that the 
quasiparticlcs and the excited states of the nucleus are defined with respect to 
the initial (final) nuclear ground state |0^ > (|0^ >). and c„ denote proton 
particle creation and neutron hole annihilation operators, respectively. 

For a complete presentation of the renormalized QRPA method we have to 
discuss also limitations of this nuclear structure method. There are no doubts 
that RQRPA offers advantages over QRPA: 

i) There is no collapse of RQRPA solution for a physical value of the particle- 
particle interaction strength. 

ii) The RQRPA Hamiltonians demonstrate better mutual correspondence like 
the QRPA ones Nevertheless the RQRPA has also several shortcomings: 
i) In the framework of the RQRPA the Ikeda sum rule (here given for 
Gamow- Teller /^-transitions) 

S^S.{l+)~S+il+) = _L_^|(^,i+|;3_|o)|2 (14) 



= 3{N-Z) (15) 

is not fulfilled. It is worthwhile to notice that the completeness of the excited 
1+ states is the only condition needed to derive the Ikeda sum rule. In the 
QRPA this sum rule is identically fulfilled, while in RQRPA the violation is 
proportional to the renormalizing matrix V, which depends on the strength of 
the residual interaction. In the framework of the RQRPA the Ikeda sum rule 
is exhausted only to 70-80%. It is supposed that the omission of the scattering 
terms in both the construction of the excited states and the evaluation of one- 
body densities could be the reason of it. At the moment there is work under 
progress including these operators in the framework of extended RQRPA |Q. 
There is a hope that in such a model the Ikeda sum rule could be restored, 
ii) The two sets of intermediate nuclear states generated from the initial and 
final ground states are not identical because of the considered QBA or RQBA 
schemes. Therefore the overlap factor of these states is introduced in the theory 
as follows fl^: 

< -^mj-^m^ >~ [<3jM:QjM ] ~ Skk'Su'^kf^t^'Ui^t.' X (16) 

fik<.i'l. fj,'k' <.v'l' 

{k,i, j)r;,i, {k,i,j)- f™; {k,i, j)r;,i, {k,i,j)), 
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Figure 1: One example diagram for a possible contribution to neutrinoless dou- 
ble beta decay in i?-parity breaking SUSY-extensions of the standard model. 



with 

w = «£4?'+4>i?'- (17) 

Here, Q"}m *3jm' ^'''^ respectively phonon annihilation and creation op- 
erators for the initial and final nuclear states. We note that in the previous 
calculations u was approximated by the unity. However, in that case the over- 
lap factor is dependent on the phases of the occupation BCS/HFB amplitudes 
u's and v's, which arc in principal arbitrary. A negligible difference between the 
results with the above overlap containing u and with the overlap without y one 
obtains only if the phases of the BCS/HFB amplitudes are chosen so that u is 
positive for each level. 



3 /?j3USY induced Oz//?/?-decay. Matrix elements 
and experimental constraints. 

i?-parity is a discrete multiplicative symmetry defined as Rp = (— i)3-B+l+2S^ 
where S, B and L are the spin, the baryon and the lepton quantum number. 
This symmetry is conserved in the minimal supersymmetric models (MSSM). 
A consequence of this symmetry would be, that SUSY-partners can only be 
produced in associated production and the lightest SUSY particle is stable. 

However, _R-parity might not be conserved. We consider presently popular 
models with the explicit i?p-violation. The i?p- violating part of the superpoten- 
tial breaking lepton number conservation and relevant to Oj^/3/3-decay is 

Win, = KjkLiQjDk. (18) 

Here L, Q arc lepton and quark doublets while U, D are lepton and up, 
down quark singlet superfields. Indices i,j,k denote generations. An example 
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Figure 2: Two-nucleon (a), one-pion (b) and two-pion (c) mode of neutrinoless 
double beta decay in i?-parity violating SUSY extensions of the standard model. 



of the diagram contributing to 0z//?/3-decay is shown in Fig. [l[ It involves the 
lepton number violating interactions originated from the superpotential (p^). 

The supersymmetric model only gives the underlying transition of a down- 
quark to an up-quark. Then this transition has to be transformed to one go- 
ing from a neutron to a proton. Looks natural and most straightforward to 
incorporate the quarks in the nucleons. In this way, we come up with the 
well known two-nucleon mode. The corresponding effective I/ip SUSY induced 
nucleon-nucleon interaction is shown in Fig. |^a). But the intermediate SUSY- 
partners are very heavy particles, and therefore in the two-nucleon mode the 
two decaying neutrons must come very close to each other, what is suppressed 
by the nucleon repulsion. 

Another possibility is to incorporate quarks undergoing the ^pSUSY tran- 
sition d -\- d — > u A~ u + 2e~ not in nucleons but in virtual pions [|. There 
are two possibilities to set up the pion contributions. Namely, only one quark- 
antiquark pair ud is placed in an intermediate pion leading to a diagram shown 
in Fig. 23) or both ud pairs are placed separately in two intermediate pions, as 
in Fig. 2z). The scale on which the pion-contribution is enhanced compared to 
the two-nucleon mode is the ratio of the nuclear form factor cut-off to the pion 
mass. As this factor favors for each nuclear decay process the pion emission, 
the two-pion mode is expected to be dominant over the two-nucleon and also 
the one-pion mode. 

The half-life for neutrinoless double beta decay regarding all the three pos- 
sibilities of hadronization of the quarks can be written in the form 



(19) 



Gor(^ 

TO 



■M 



2N 



Here Gqi is the standard phase space factor tabulated for various nuclei in 
Ref. jnl and = 850 MeV. The factors -qq and 77^ are ruled by the SUSY- 



8 



parameters and by the strength of i?-parity violating couphng. This factors will 
be constrained by the non-observation of neutrinoless double beta decay, and 
under reasonable assumptions for the SUSY-parameters (especially the SUSY 
masses), a limit on the i?-parity violating couphng X[^^ can be derived. The 
nuclear matrix elements -MV^j governing the two-nucleon mode were presented in 



Ref. 



. The one- and two-pion modes contribute to the nuclear matrix element 
which consists of the four partial matrix elements 
m 



M 



ttN 



[MgT,!^ + Mt.Itt) + a^'" (AfGT,27r + Mt,2 



(20) 



two (Gamow- Teller and Tensor type) elements for each pionic mode. The two 
types of matrix elements are given by the expression 



M, 



GT,k7T — 



Mt 



kn 



with fc = 1, 2 



(21) 



(Ti ■ (Jo 



with = rriTrrij, Vij = \fi — fj\ and — {fi — The structure functions 

Fi''\xTr) and Fj'^'' (x-n-) have their origin in the integration of the propagator of 
the intermediate particles and take the following form: 



f[^\x) = [x-2)t 



F^^\x) = [i + ix + x'f— 



F^^>{x) ^ {x + l)e~ 



(22) 
(23) 



The structure coefficients a-^^ and a'^^ depend on the hadronization of the 
quarks. For the two-pion coefficient a^^ a hadronic matrix element of the type 
(tt"*"! JjJi 1 7r~) {i — P,S,T) for the hadronic currents Ji needs to be evaluated. 
This can be done in either a non-relativistic quark model (QM) or by inserting 
the vacuum state in between the two currents (vacuum insertion approxima- 
tion, VIA). As the VIA neglects contributions from other intermediate states 
than the vacuum, it will give a more conservative limit. The values for the 
structure coefficients a'^'' = -4.4 • 10"^ ^nd a^'^ 0.2(VIA), 0.64(QM) reveal 
the above made statement on the dominance of the two-pion over the one-pion 
mode. Further it will be seen that not only the structure coefficient but also 
the corresponding nuclear matrix elements MoT^k-n and Mt^r-k of the two-pion 
mode is larger than those of the one-pion mode. 

For further calculation the nuclear matrix element governing the SUSY- 
mechanism of 0!//3/3-decay are transformed to ones containing two-body matrix 
element in relative coordinate. One arrives at the expression: 



< Oi2 > = 



E 

klkl 



(-)■" 



{2J+1) 



3k 
ji' 



31 

jk' 



J 
J 
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Table 1: Nuclear matrix elements for SUSY two-nucleon, one-pion and two-pion 
mechanisms of neutrinoless double beta decay for ^^G'e(0"'') — Se{0~^) nuclear 
transition within the pn-RQRPA and fuU-RQRPA. The 12-level (the full 2-4hio 
major oscillator shells) model space has been considered. The presented values 
has been calculated for gpp = 1.0. The BM, QM and VIA denote bag model, 
non-relativistic quark model and vacuum insertion approximation, respectively. 

mechan. two-nucleon 



M.E. 


McT,N 


Mf,n 


Mgt' 


Mf' 


Mt' 


Mq 






[10-2] 


[10-2] 


[10-2] 


[10-3] 


[10-3] 


(BM) 


(BM) 














(NQM) 


(NQM) 


pn- 


7.05 


-2.48 


-1.04 


3.76 


-2.38 


-116 


5.7 


RQRPA 












-156 


2.7 


full- 


4.32 


-1.26 


-0.65 


1.94 


-0.76 


-75 


5.5 


RQRPA 












-95 


4.2 


mechan. 




one-pion 






two-pion 




pion 


M.E. 


MgT,Itv 


Mf.m 


Ml. 


MgT,27t 


Mf.2^ 


M2. 
















(VIA) 


(VIA) 














(QM) 


(QM) 


pn- 


1.296 


-1.023 


-22 


-1.341 


-0.653 


-1.99 


-754 


RQRPA 














-2364 


fuU- 


0.840 


-0.450 


-31 


-0.810 


-0.385 


-439 


-470 


RQRPA 












-1403 


-1434 



< Oj^ II [Cpfc'5ni']j II >< rmf\rm, >< .rm, II [c+ c„,]j || 0+ > 

X <pk,pk'-J\f{ri2)TtT+Oi2f{ri2)\nlM'\J> ■ (24) 
with a short-range correlation function 

/(r) ^ 1 (1 -6r2) with a==l.lfm2 and 6 = 0.68fm2, (25) 

which takes into account the short range repulsion of the nucleons. 

The calculated nuclear matrix elements for the Qv(3(3-Aecsy of A=76 isotope 
within the pn-RQRPA and the fuU-RQRPA are presented in Table 1. The con- 
sidered single-particle model space has been the 12-level model space (the full 
2 — Afiui major oscillator shells) introduced in Ref. [^. One should note that 
the calculation of the matrix elements in QRPA collapses with increasing in- 
teraction strength, while a calculation with RQRPA stays stable in the whole 
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range of physical interaction strength. The nuclear matrix elements listed in 
Table 1 have been obtained for the gpp = 1 {gpp- parameter used to renormal- 
ized particle-particle interaction of the nuclear Hamiltonian) . By glancing the 
Table 1 we note less important role of the one-pion exchange SUSY mechanism 
and that the two-pion exchange nuclear matrix elements clearly dominate over 
nuclear matrix elements of the two-nucleon mechanism. There are two sources 
of this enhancement. The first source has pure nuclear origin as the poten- 
tial determined by the pion-exchange with mass about 140 MeV is favored in 
comparison with the potential determined by the cut-off with value about 850 
MeV. The second source of enhancement has its origin in the hadronization 
of the ]jtp SUSY effective vertex operator u'-f^d ■ u'-f^d ■ ePae'^ replaced by its 
hadronic image tt^ • ePfje'^. The enhancement occurs due to the coincidence of 
the pseudoscalar quark bilinears u'^^d with 7r-meson field. 

For deducing constraints on the i?-parity violating SUSY parameters from 
the non-observation of Oi^/9/3-decay we shall use values of nuclear matrix elements 
obtained within the fuU-RQRPA. The current experimental lower bound on the 
^^Ge Or//3/3-decay half-life (l| is 

^Ow3/3-exp^Q+ ^ > 1.1 X 1025 years 90% c.l. (26) 

With the above calculated nuclear matrix elements lower limit can be trans- 
formed to the following upper limits for the 1st generation Yukawa coupling 
constant 

^ '■^<°-«''»'*(To5^)'(Tosfc7)"' m 

We point out at that the uncertainties of the nuclear structure calculations are 
smaller than those from the hadronic matrix elements. The limits in Eqs. (|27|)- 
( ^8|) correspond to the VIA (QM) calculations of the hadronic matrix element. 

These limits are much stronger than those previously known and lie beyond 
the reach of the near future accelerator experiments (though, accelerator exper- 
iments are potentially sensitive to other couplings than X'm). To constrain the 
size of X'lii one needs to make assumptions on the masses of the SUSY-partners. 
If the masses of the SUSY-partners would be at their present limit ||l6| , one could 
constrain the coupling to X'^^^ < 6.0(3.3) • 10"^. A conservative bound can be 
set by assuming all the SUSY-masses being at the " SUSY-naturalness" bound 
of 1 TeV, leading to A'^ < 8.2(4.4) • lO^^. 

These results show that the non-observation of 0i/(3f3-deca.y strongly limits 
extensions of the standard model of electroweak interaction. Although a many- 
body problem needs to be solved the improvement of the limits is so large that 
it overcomes the uncertainties in the nuclear and hadronic matrix elements and 
leads to limits that are much stronger than those from accelerator and non- 
accelerator experiments. 
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